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ABSTRACT 

We report spectral, imaging, and variability results from four new XMM-Newton observations and 
two new Chandra observations of high-redshift (z > 4) radio- loud quasars (RLQs). Our targets 
span lower, and more representative, values of radio loudness than those of past samples of high- 
redshift RLQs studied in the X-ray regime. Our spectral analyses show power-law X-ray continua 
with a mean photon index, (r) = 1.74 ± 0.11, that is consistent with measurements of lower redshift 
RLQs. These continua are likely dominated by jet-linked X-ray emission, and they follow the expected 
anti-correlation between photon index and radio loudness. We find no evidence of iron Ka emission 
lines or Compton-reflection continua. Our data also constrain intrinsic X-ray absorption in these 
RLQs. We find evidence for significant absorption (JVh ~ 1.7 x 10 22 cm^ 2 ) in one RLQ of our sample 
(SDSS J0011+1446); the incidence of X-ray absorption in our sample appears plausibly consistent with 
that for high-redshift RLQs that have higher values of radio loudness. In the Chandra observation of 
PMN J2219— 2719 we detect apparent extended (~ 14 kpc) X-ray emission that is most likely due to 
a jet; the X-ray luminosity of this putative jet is * 2% that of the core. The analysis of a 4.9 GHz 
VLA image of PMN J2219— 2719 reveals a structure that matches the X-ray extension found in this 
source. We also find evidence for long-term (450-460 days) X-ray variability by 80-100% in two of 
our targets. 

Subject headings: cosmology: observations — X-rays: galaxies — galaxies: active — quasars 



1. INTRODUCTION 

Active Galactic Nuclei (AGN) are tracers of super- 
massive black hole (SMBH) growth in the centers of 
galaxies. This activity can be used to find and probe 
the evolution of massive galaxies at cos mological epoch s 
close to the period of reionization (e.g.. lFan et al.ll2006h . 
The study of high-redshift quasars (z > 4) is particu- 
larly important because it provides information on the 
period when the first SMBHs were formed in the Uni- 
verse. X-ray studies reveal conditions in the immedi- 
ate vicinity of these SMBHs as well as their larger scale 
environment (e.g., with intrinsic or intervening absorp- 
tion). Most X-ray studies of high-redshift quasars have 
focused on radio-quiet quasars (RQQs; those with radio- 
loudness 10 parameter Rl < 10). The main result of 
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these studies is that high-redshift RQQs have similar 
X-ray properties (e.g., X-ray power-law photon index, 
r, and optical-to- X-ray power-law slope 11 , a ox ) to those 
at lower redshift (z = 1 — 3) of matched luminosity (e.g., 
IShemmer et "all [20061 : Dust et al.1 12007I) . 12 Thus, high- 
redshift RQQs seem to be feeding and growing in a gen- 
erally similar manner to those at lower redshift. 

Radio-loud quasars (RLQs; those with Rl > 10 — 100) 
have strong radio emission which is usually linked to the 
presence of jets emitted from the AGN at relativistic 
speeds (e.g.. iBegelman et al.lll984l) . X-ray spectral and 
variability studies of high-redshift RLQs provide a probe 
of their pc-scale jet-linked X-ray emission, since this 
jet-linked emission often dominates the observed X-ray 
luminosity. Due to this jet-linked component, RLQs 
have stronger X-ray emissi on than RQQs of simi lar op- 
tical/UV luminosities (e.g.. IZamorani et al.l ll981). This 
exces s in X-ray emission in creases with the radio loudness 
(e.g., IWorrall et al.l lL987h : for example, for log Rl ~ 4 
the X-ray luminosity of RLQs can be up to ~ 10 times 
higher than the luminosi ty of RQQs with s imilar opti- 
cal/UV luminosity (e.g., I Miller et all 1201 1| ). The jet- 
linked emission of RLQs also ma kes their X-ray spec- 
tra harder as Rl incre ases (e.g., iWilkes fc Elvisl U987t 
iReeves fc Turner! 120001 ) . RLQs re present 10-20% of 
the whole quasa r population (e.g., Ilvezic et"aT1 l2002t 
Uiang et"aT]|2007| ). Notably, the frequency of radio jets in 
luminous quasars seems to change significantly with red- 

11 The optical-to- X-ray power-la w slope is defined as q o x = 
log(/ 2 kcv//2500A)/log(^2kev/j y 2500A) l ITananbaum et al.l ll979h. 

12 For RQQs, a ox sho ws a dependence upon the UV luminosity 
(e.g., ISteffen et al.N2006ll . For these sources ther e is also evidence 
that r may depend on the X-ray luminosity (e.g. . [Saez et al.l 2008) 
and/or the accretion rate (in terms of the Eddington ratio; e.g., 
IShemmer et aI.ll2008T) . 



2 



shift and UV luminosity fe.g. JJiang et al.ll2007[ ). Thus, it 
is of interest to assess if the physical processes operating 
within jets also change with redshift; any such changes 
should probably be most apparent at the highest red- 
shifts. RLQs at z > 4 have been studied much less well 
in the X-ray regime than correspondingly high-redshift 
RQQs. The number of sensitive exploratory X-ray ob- 
servations (and thus X-ray detections) is much smaller 
(20 — 30 in total), and the number of RLQs with quality 
X-ray spectra is even smaller still (< 10). Those z > 4 
RLQs that have been targeted for X-ray sp ectroscopy are 
usually highly radio lou d (Rl > 1000; e.g. JWorslev et all 
120041 lYuan et al.ll2006[ ). Such highly radio-loud objects 
make up only a small fraction (< 5%) of the RLQ pop- 
ulation, so they are not representative objects. 

X-ray observations of high-redshift RLQs with good 
angular resolution can also facilitate study of kpc-scale 
X-ray jets. Two of the most popu lar models for this kpc- 
scale emission (e.g., IWorrallll2009l and references therein) 
are (1) an extra high-energy synchrotron component, and 
(2) inverse Compton (IC) scattering by relativistic elec- 
trons of photons of the cosmic microwave background 
(CMB; or the IC/CMB model). For synchrotron mod- 
els, we do not obviously expect that the X-ray emission 
properties would depend on redshift. The energy den- 
sity of the CMB increases as (1 + z) 4 ; as a consequence, 
an implication of the IC/CMB model is that for high- 
redshift RLQs (z > 4) the X-ray emission from extended 
jets may outshine that from the core (e.g., ISchwartzl 
I2002T) . High-redshift sources with such dominant X-ray 
emiss ion from extended jets have not yet been reported 
(e.g.. lLopez et al.l 12006). This IC/CMB enhancement 
of the X-ray emission could also be observable via no- 
ticeable systematic changes with redshift of a x and the 
radio-to-X-ray power-law slo pe 13 , a rx \ such increments 
have not been obs erved (e.g.. iMarshall et al.ll2005l 120111 ; 
iMiller et al.ll20Ll . 

The potential X-ray absorption in representative high- 
redshift RLQs is also a critical topic to study. Changes 
in the strength and frequency of RLQ intrinsic X-ray ab- 
sorption (i.e., TVh) with redshift have been discussed by 
many authors (e.g . , lElvis et al.lll994at [Cappi et al1ll997t 
Fiore et al.l 119981 : iReeves fc Turner! l2000t iPage et all 



spect ra show evidence for substantial X-ray absorption 
/e.g.- lFabian et al.ll200lUWorslev et al.ll2004HPage et al.l 



l2005i Yuan et al.l |2006|) . a finding that is broadly con- 



20051 : lYuan et al.ll2006f ) . These studies demonstrate that 
the fraction of RLQs showing low-energy X-ray cutoffs 
appears to rise with redshift. These cutoffs are typically 
attributed to the presence of intrinsic X-ray absorption. 
Column densities of ATh ~ (1-3) xlO 22 cm" 2 are often 
seen at z > 2. RQQs do not show clear evidence for 
increasing X-ray absorption with redshift; this contrast 
with RLQs demonstrates that there must be a connec- 
tion between X-ray absorption and radio loudness. The 
absorbing gas is thought to be associated with the envi- 
ronments of RLQs, but its precise nature is unclear: it 
may be circumnuclear, located in the young host galaxy, 
or entrained by the radio jets. Based on Swift/GRB 
data set it has also been proposed recently that this ab- 
sorption m ay arise in the di ffuse intergalactic medium 
(IGM; e.g.. lBehar et al.ll201lD . Three of the seven highly 
radio-loud (log Rl > 3) quasars at z > 4 with X-ray 

13 The radio-to-X-ray power -law slope is d e fined as a rx = 
l°g(/ikcv//r)/log(i'i] J0 v/i y r); in Marshall ct al. (2011) the radio 
frequency used is v r = 8.64 GHz. 



sistent with an extrapolation of the trend observed at 
lower redshifts. However, the source statistics are lim- 
ited, and it is precarious to generalize the results from 
a few extreme highly radio-loud objects to the high- 
redshift RLQ population as a whole. Note also that 
moderate low-energy absorption could be mimicked by 
a curved intrinsic spectr um; this possibility c annot be 
ruled out for RLQs (e.g.. iTavecchio et al.ll2007D . 

In past work, we have targeted representative flat- 
spectrum RLQs at z fa 3.5 — 5, typically having moder- 
ate Rl values (Rl ~ 9 — 400), using Chandra snapshot 
obser vations (5-10 ks; iBassett et al.l 120041 : lLopez et al.l 
|2006| ). These snapshot observations gave information 
about the level of X-ray emission and allowed basic pop- 
ulation studies. However, they did not allow detailed 
physical studies of individual objects. Thus, in order 
to investigate some of the issues described above, we 
have targeted selected objects from our snapshot sam- 
ples for longer follow-up observations with XMM-Newton 
and Chandra. These observations provide spectral con- 
straints as well as information about long-term X-ray 
variability. Furthermore, the observations with Chan- 
dra allow effective searches for kpc-scale X-ray jet emis- 
sion. We were awarded time to observe four targets with 
XMM-Newton and two targets with Chandra. The two 
targets observed with Chandra were selected to follow- 
up hints of X-ray extension based upon analyses of the 
snapshot data. We supplemented the follow-up Chandra 
imaging with high-resolution Very Large Array observa- 
tions to constrain the presence and extent of putative 
jets. 

The layout of this paper is as follows: in §2 we de- 
scribe our sample, the X-ray data reduction, and the ra- 
dio data reduction; in §3 we provide a discussion of the 
most important results found in this work; and in §4 we 
summarize our results and main conclusions. Through- 
out this paper, unless stated otherwise, we use cgs units, 
the errors listed are at the la level, and we adopt a 
flat A-dominated universe with Hq = 70 km s _1 Mpc -1 , 
Q A = 0.7, and Q M = 0.3. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. X-ray data 

In this work we perform analyses of four new 
XMM-Newton observations and two new Chandra obser- 
vations of high-redshift (z = 3.6 — 5.0) RLQs with moder- 
ate Rl- These observations were obt ained as a follow-up 
of Chandra snapshot observations (IBassett et all 12004 
lLopez et al.l 120061 : IShemmer et~aTl |2006j). The original 
snapshot observations were also reprocessed with the 
main objective of studying any long-term X-ray vari- 
ability. The observation logs of the sources analyzed 
in this work, that include observation dates, observed 
total counts, total exposure times, and observational 
identification numbers, are presented in Table [1] for the 
XMM-Newton observations and Table [2] for the Chan- 
dra observations. Additionally, in Table [3] we present 
basic optical, X-ray, and radio properties for our sample; 
among these are the redshift and Rl of each source. The 
radio-loudness parameters range from w 90 — 7000 with 
five out of the six sources having Rl < 200. From Ta- 
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TABLE 1 

Log of XMM-Newton Observations 











EXPOSURE 


TIME b (ks) / TOTAL COUNTS 


OBJECT NAME 


Q 2000.0 a 


^ooo. o a 


OBS. DATE / ID 


MOS1 


MOS2 


pn 


SDSS J0011+1446 
PMN J0214-0518 
SDSS J0839+5112 
SDSS J1309+5733 


00 11 15.23 
02 14 29.30 
08 39 46.21 
13 09 40.69 


+ 14 46 01.8 
-05 17 44.6 
+51 12 02.8 
+57 33 09.9 


2010Jan07/0600090101 
2010Jan01/0600090401 
2006Aprl2/0301340101 
2006Apr22 /0301340501 


29.3(38.5)/761±28 
58.9(60.5)/186t^ 
9.6(26.6)/lQl±iJ 
29.4(54.6)/85tg° 


29.3(38.5)/773t28 
54.3(60. 5)/205ti4 

9.6(26.6)/87lg° 
28.6(54.6)/112tJj 


24.5(36. 4)/2336t4g 

46.0(58.9)/773±28 
7.4(24.9)/314iig 

18.3(52.9)7235^^ 



Note. — Errors on the X-ray counts were computed according to Tables 1 and 2 of Gchrels (1986). 
a Optical p ositions in J2000.0 equatorial co ordinates. The positions of SDSS J0011+1446 and SDSS J0839+5112 and SDSS J1309+5733 are 
presented in Adclman-McC arthv et al.l J2007fl ■ and the position of PMN J0214-0518 is presented in lMcMahon et al.l <2002fl . 

The exposure times and photon counts are obtained after screening the data for flaring. The exposure times before the screening are in parentheses. 
c Total source counts in a circular region of radius of 30" centered on the X-ray source. For the MOS cameras and pn camera the photon energies 
in consideration arc the 0.5—10 keV and the 0.3—10 keV observed band respectively. 



TABLE 2 
Log of Chandra Observations 



OBJECT NAME 



02000.0 



EXPOSURE TIME 
OBSERVATION DATE/lD (ks) 



TOTAL COUNTS" 



Ref. d 



New observations 



PMN J0235- 
PMN J2219- 



1805 
2719 



02 35 02.50 
22 19 35.32 



-18 05 51.0 
-27 19 03.3 



2009 Feb 22 / 10306 
2009 Sep 12 / 10305 



19.9(21.5) 
37.9(38.5) 



141 

1803 



12 
+43 



Archival observations 



SDSS J0011+1446 


00 11 15.23 


+14 46 01.8 


2003 May 20 / 3957 


3.5(3.8) 


135±g 


2 


PMN J0214-0518 


02 14 29.30 


-05 17 44.6 


2003 Nov 26 / 4767 


4.0(4.3) 


44 ±y 


3 


PMN J0235-1805 


02 35 02.50 


-18 05 51.0 


2004 Feb 07 / 4766 


3.9(3.9) 


20+5 


3 


SDSS J0839+5112 


08 39 46.21 


+51 12 02.8 


2003 Jan 23 / 3562 


4.9(5.0) 


82+1° 


4 


SDSS J1309+5733 


13 09 40.69 


+57 33 09.9 


2003 Jul 12 / 3564 


4.6(4.7) 


311? 


i 


PMN J2219-2719 


22 19 35.32 


-27 19 03.3 


2003 Nov 19 / 4769 


8.1(8.2) 


2261 Jl 


3 



Note. — Errors on the X-ray counts were computed according to Tables 1 and 2 of [Gehrcls ( 198^). 
a Optical positions in J2 000.0 equatorial coordinates. The positions of PMN J0235-1805 and PMN J2219-2719 are presented 
in lMcMahon et al.1 1 120021) . 

The exposure times and photon counts are obtained after the screening of the data. The exposure times before the screening 
are in parenthesis. 

c Total source counts with energies in the 0.5-8 keV observed band in a circular region of radius of 3.9 arcscc (8 pixels) centered 
on the X-ray source. 

d References: (1) This work; (2l lShemmer et al.1 J2006I) ; O HLopez et al.l d2006f) : C4) IBassett et al.l 120041) . 



ble[3]we also notice that our RLQ sample is on average 
~ 6 times brighter in X-rays tha n RQQs o f similar UV 
luminosity (based on Uust et all 2007 and iMiller et al.l 
2011). This is suggestive evidence that the X-ray con- 
tinua of our sample are dominated by jet-linked emis- 
sion. In Figure Q] we have plotted our sources (with 
stars) in a diagram of log Rl vs. z. We also have 
adde d to this plot RL Qs from other high-redshift sam- 
ples (page et al.1 120051; IShemmer et al.ll2005l [Yuan et al.1 
2006: lYoung et al.l l2009f ) with high-quality X -ray spec- 



tra. Our sample is concentrated toward high redshifts 

14 In the added sample of RLQs there are some sources that have 
been analyzed in more than one work. The sources with overlap are 
PMN J0525-3343, RX J1028. 6-0844, SPSS J1430+4204 (these 
three sources were analyzed in lPage et al.|[ 2005; Yuan et al.H2006l) 
and S PSS J1510+5702 (analyzed in IPage l et al 1120051 ; lYuan et all 
2006; lYoung et al.l [20091 ). We avoid repeated sources m our ex- 
tended sample by s electing first l y the spectral properties of the 
so urces presen t ed in IPage et al.l 120051) , secondly th ose presented 
inlYuan et al.l 020061 ). and thirdly those presented in Young et al.l 
H2009I) . 



and typically representative radio-loudness parameters, 
where it significantly improves the source statistics. This 
can be seen from Figure [1] by comparing the R l of our 
sources to the Rl distribution fo und from SDSS/FIRST 
detections bv llvezic et al.l |2002) (dotted-line curve). 

The XMM-Newton data were analyzed with the Sci- 
ence Analysis Software (SAS) version 10.0.0 provided 
by the XMM-Newton Science Operations Centre (SOC). 
The event files were obtained by including the events 
with flag = 0, pattern < 12 (pattern < 4), and 200 < 
PI < 12,000 (150 < PI < 15,000) for the MOS (pn) 
detectors. The event files were also temporally filtered 
in order to remove periods of flaring activity. 15 This 
flaring activity tends to saturate the instrumental back- 
ground and therefore reduces the good-time intervals. 

15 To filter the event files of periods of flaring activity, we create 
light curves using single events with energy greater than 10 keV. 
Using these light curves we find and remove from the event files 
the time periods where the count rate is above ~0.5 count s _1 for 
the MOS cameras and ~1 count s —1 for the pn camera. 
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TABLE 3 

Optical, X-ray and Radio Properties of Surveyed Radio-Loud Quasars at High Redshift 



OBJECT NAME 


z 


ABl450 a 


Sl.4GHz 

(mJy) 




aox d 


log R L C 


log L 5 GHz f 


log Z/2500 S 


log i2kcV h 


Alog Z/2kcv' 


SDSS J0011+1446 


4.97 


18.1 


36 


-0.12 


-1.31+0.03 


2.01 


34.21 


32.08 


28.66 


1.10 


PMN J0214-0518 


3.99 


18.8 


41 


+0.65 


-1.53+0.04 


2.25 


34.01 


31.65 


27.67 


0.41 


PMN J0235-1805 


4.31 


18.9 


36 


+0.26 


-1.50+0.04 


2.28 


34.05 


31.64 


27.74 


0.49 


SDSS J0839+5112 


4.39 


18.8 


43 


+0.14 


-1.40+0.03 


2.33 


34.16 


31.71 


28.07 


0.77 


SDSS J1309+5733 


4.27 


19.3 


11 




-1.51+0.03 


1.98 


33.58 


31.48 


27.53 


0.39 


PMN J2219-2719 


3.63 


20.3 


304 


-0.26 


-1.09+0.03 


3.85 


34.93 


30.96 


28.11 


1.34 



a The monochromatic AB magnitude at r est-frame wavelength 1450 A. For PMN J0214-051S, PMN J0235-1805, and PMN J2219-2719 we used the 
R magnitudes found in lHook et al.1 i2002fl to calculate AB1450 by using the empirical relation AB1450 — R — 0.684z + 3.10. For SDSS J0011+1446, 
SDSS J0839+5112, and SDSS J1309+5733 we used the i-band SDSS magnitudes to calculate AB145Q by using the empirical relation AB1450 — i — 0.2. 
Notice that these estimates of AB\4<$q are reliable within £3 0.1—0.2 mag for the redshifts of our sources. 

b The 1.4 GHz fl ux density; from the NRAO VLA Sky Survey (NVSS) JCondon et all [19981) . For sources covered by VLA FIRST Sky survey 
(White et al. 1997), their 1.4 GHz flux density differ in less than 10% with the values presented. 

c The radio power-law slope between 1.4 and 5 GHz (observed frame). The 5 GHz flux density is from the Green Bank 6 cm (GB6; Gregory ct al. 
1996) or Par kes- MIT- NRAO (PMN; Griffith & Wright 1993) surveys. 



The optical-to-X-ray power-law slope a D 



log(/2kev//25ooA)/l°s( t/ 2kcv/^250oA) ' The average difference between measured and predicted a ox for 



RQQ is (Aq ox ) 0.3, based on lSteffen et aU 112006ft ^ox-^sooA relation. 

e Radio-loudncss parameter, defined as Rl — fsGHz //4400A ( rcs t frame; Kcllcrmann ct al. 1989). The flux densities at 5 GHz in the rest-frame were 
obtained from the flux densities at 1.4 GHz in the observed-frame assuming a power-law spectrum with spectral slope given by a r ; when a r is not 
available we assume a. r — 0. The fluxes at 4400 A in the rest-frame are obtained from the AB magnitudes assuming a power-law spectrum with 
spectral slope of —0.5 (Vandcn Berk ct al. 2001). 

^ Logarithm of the monochromatic luminosity at 5 GHz. The luminosities were obtained from the flux densities at 1.4 GHz in the observed-frame 
assuming a power-law spectrum with spectral slope given by a r ; when a r is not available we assume a r — 0. 

^ Logarithm of the mono chromatic luminosity at 2500 A. These were computed from the A.B1450 assuming and optical power-law of a — —0.5 
dVanden Berk et al.|[200ll) . 

Logarithm of the absorption-corrected monochromatic luminosity at 2 kcV. 
1 The difference between measured and the predicted log L,2kcV f° r RQQs; based on the Uust et al relation 
log i 2 kcv = 0-709 x log L2500 + 4.822 . 
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Fig. 1. — Radio-loudness parameter versus redshift for our 
sample (stars) and z > 2 RL Qs observed with Chandra and 
XMM-Newton (circles are from Pa ge et al.l 20 05, triangles from 
IShemmer eTal 2005. squar es from lYuan et aLl 2006. and inverted 



triangles from lYoune et aLl 2009). The dotted curve along the left 
side of the p l ot sho ws the relative number of RLQs versus Rl from 
llvezic et aLl l|2002l '). Note that our sample significantly improves 
the source statistics for moderately radio-loud quasars at the high- 
est redshifts. 

The flaring is especially strong in the observations of 
SDSS J0839+5112 and SDSS J1309+5733, where the ef- 
fective observation time has been reduced by 40% and 
50%, respectively, when compared with the original ex- 
posure time (see Table [lj . For the three EPIC detec- 
tors (MOS1, MOS2, and pn), the source and background 



spectra of each quasar were extracted from a circular re- 
gion with an aperture radius of 30" 16 and an annular 
source-free region with a 50" inner radius and 100" outer 
radius, respectively. The numbers of counts in the source 
regions for the XMM-Newton cameras are presented in 
Table [T] These counts are selected at energies of 0.5- 
10 keV for the MOS cameras and 0.3-10 keV for the pn 
camera. The combined numbers of counts in the EPIC 
cameras range approximately from 500 to 4000 counts. 
The response hies used to analyze the source spectra were 
created with the SAS tasks rmfgen for the redistribu- 
tion matrix files (RMFs) and ARFGEN for the ancillary 
response files (ARFs). 

The Chandra observations (new and archival) were an- 
alyzed using the standard software CI AO 4.2 provided by 
the Chandra X-ray Center (CXC). The telemetry format 
used was Faint mode for all the observations with the 
exception of obsID 10306 for which Very Faint mode 17 
was used. Standard CXC threads were employed to 
screen the data for status, grade (ASCA grade 0, 2, 3, 
4 and 6 events), and time intervals of acceptable as- 
pect solution and background levels, only. The source 
and background spectra of each quasar were extracted 
(using the CIAO task dmextract) from a circular re- 
gion with an aperture radius of 3.9" 18 (8 pixels) and an 
annular source-free region with an inner radius of 4.9" 
(10 pixels) and an outer radius of 39" (80 pixels), re- 

16 The 1.5 keV encircled energy at 30" from the center of a 
point source is between 80-90% for the EPIC cameras. Note that 
throughout this work we have applied an aperture correction factor 
when we calculate fluxes or luminosities. 

17 Very Faint mode offers the advantage of reduced background 
after ground processing. 

18 The 1.5 keV encircled energy at 3.9" from the center of a 
point source is ~ 99% for ACIS-S . 
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Fig. 2.— Full band (0.5-8 keV) Chandra images of PMN J0235-1805 (left panel) and PMN J2219-2719 (right panel). Each panel 
spans 4.4" X 4.4" on the sky; North is up, and East is to the left. Additionally, in each panel the VLA contour plots at 4.9 GHz are shown. 
The image resolutions of the radio contour plots are 0.55" X 0.30" at PA = -7.7° and 0.80" X 0.32" at PA = 7.5° for PMN J0235-1805 
and PMN J2219— 2719, respectively. The contour plots begin at 0.4 mjy per beam and increase by factors of 4. Notice that in each panel 
there is an offset between the X-ray image and the radio contours of about 0.2" — 0.3"; this offset is within the expected astromctric error 
of Chandra. The elliptical restoring beams are plotted in the lower-right inset of each panel. 



spectively. The numbers of counts in the source regions 
at energies of 0.5-8 keV are presented in Table [5] The 
0.5-8 keV source counts in the new Chandra observations 
of PMN J0235-1805 and PMN J2219-2719 are approx- 
imately 140 and 1800 counts, respectively (see Table [2]). 
The 0.5-8 keV source counts in the archival observations 
range from w 20 — 230 counts (see Table[5|). The response 
files used to analyze the source spectra were created with 
the CIAO tasks MKACISRMF for the RMFs and mkarf 
for the ARFs. 

2.2. Radio data for PMN J0235-1805 and 
PMN J2219-2719 

The new Chandra-observed targets, PMN J0235-1805 
and PMN J2219-2719, were observed at the NRAO 
Very Large Array (VLA) in A configuration, for a nom- 
inal resolution of 1.0" and 0.39" at 1.4 and 4.9 GHz, 
respectively. PMN J0235— 1805 was observed in Aug 
2007 (observation proposal code AL696), with 100 min- 
utes spent at 1.4 GHz and 28 minutes at 4.9 GHz. 
PMN J2219-2719 was observed on 16 Oct 2004 (ob- 
servation proposal code AC755), for 35 minutes at 1.4 
GHz and 12 minutes at 4.9 GHz. Observations were ob- 
tained in VLA continuum mode, with two intermediate- 
frequencies (IFs) and full polarization, for a total band- 
with of 86 MHz. The data were calibrated using the ab- 
solute flux-density calibrators 3C48 and 3C147 and the 
secondary calibrators 0240-231 and 2248-325. Data 
editing and calibration were carried out in the Astro- 
nomical Image Processing System (AIPS). The quasars 
were imaged using uniform weighting for maximal reso- 
lution, and were self-calibrated with a phase-only solu- 
tion interval of 10 sec and amplitude solution intervals 
of 5 min. The VLA synthesized beam at 1.4 GHz and 
4.9 GHz at the position of PMN J0235-1805 has dimen- 
sions of 2.04"xl.03" at PA = -20.2° and 0.55"x0.30" 
at PA = —7.7°, respectively. The VLA synthesized 
beam at 1.4 GHz and 4.9 GHz at the position of 



PMN J2219-2719 has dimensions of 2.70" x 1.00" at 
PA = 3.8° and 0.80"x0.32" at PA = 7.5°, respectively. 

3. RESULTS AND DISCUSSION 

Here we describe and discuss the main results obtained 
in this work. We start in SI3.1I by assessing evidence for 
extended jets in the Chandra images. We continue in i)3.2l 
by providing the main results from our spectral analyses 
with XMM-Newton and Chandra. In ^3.31 we constrain 
long and short-term variability. Finally, in SI3.4I we per- 
form correlation analyses between several of the param- 
eters that describe our spectral data; in particular we 
focus on correlations found in low redshift samples, as 
e.g., r with log Rl and Ah with z. 

3.1. Imaging analysis with Chandra 

The outstanding angular resolution of Chandra 
(~ 0.5") allows us to study any jet that has an X-ray 
spatial extension larger than ~ 5 kpc in the RLQ rest 
frame. The analysis of the snapshot obs ervations of 
PMN J0235- 1805 and PMN J2219-2719 bv lLopez et all 
( 2006) hinted that these two sources might present ex- 
tended jets in X- rays. The putative jets found by 
iLopez et al.l ([2006) were extended (1—2") in the west- 
ern and south-western direction for PMN J0235— 1805 
and PMN J2219-2719, respectively. In both sources the 
X-ray luminosity of the e xtended emission w as found to 
be « 3% that of the core (jLopez et al.ll2006f ) . 

In Figure [5] we present the new full band X-ray 
(0.5-8.0 keV) images of PMN J0235-1805 and 
PMN J2219-2719. In this figure we also have overlaid 
the contour plots of their 4.9 GHz emission based on the 
VLA data. From a simple inspection of these images 
we do not find any obvious evidence for highly extended 
X-ray emission in these sources. To check further for any 
evidence of extended emission we have compared a radial 
profile extracted from each of the Chandra images using 
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Fig. 3. — 0.5—2 keV (left panel) and 2—8 keV (right panel) radial profiles of the PMN J2219— 2719 core region, constructed from circular 
sectors shown in the upper-right inset of each panel. The radial profiles are obtained assuming positive distance in the southwestern 
direction. For presentation purposes, in the upper-right inset of each panel, we have masked the photons coming from a one arc-second 
radius circular region centered at the centroid of the source. The dashed curve shows a scaled PSF generated with MARX with energies 
close to the average energy of the photons encircled at 2" from the centroid of the X-ray image. 



circular annuli with a point-spread function (PSF) cal- 
culated with MARX. 19 The images for each source were 
made using photons in the soft 0.5-2 keV band and the 
hard 2-8 keV band. These images were obtained by 
removing pi xel randomizati on and applying the subpixel 
algorithm of lLi et ail (|2004) . For each source we obtain a 
MARX simulation assuming a point source with an X-ray 
spectrum given by the best-fitted absorbed power-law 
model from the new observations (see EI3.2I) . 

From the PSF analysis with MARX we found a hint 
of spatial extension in PMN J2219-2719 along the 
south-west direction. In Figure |3] the soft band (left 
panel) and hard band (right panel) radial profiles of 
PMN J2219-2719 along the south-west direction are 
presented. The possible X-ray spatial extension is seen 
in the soft band 20 through a count overdensity in the 
data with respect to the PSF around 2" from the cen- 
troid of the source. This count overdensity is seen in 
two data points with significantly higher counts (~ 2a) 
than those predicted from the PSF (see left panel of Fig- 
ure |3|) . To obtain a basic estimate of the level of pu- 
tative jet emission in PMN J2219— 2719 we measure, in 
the full band (0.5-8 keV), the excess of counts, with re- 
spect to a point source, that lie at >1" from the core. 
For this purpose the radial profile of the source is com- 
pared with a MARX simulation of a point source. Our 
results indicate that we have 381^ (errors from Gehrels 
1986) additional counts 21 in the data when compared 
with the simulated PSF over a circular aperture ex- 

19 MARX is a Chandra ray-tracing simulator (for details see 
|http: / /space.mit.edu/CXC /M ARX) . 

The soft band rest-frame coverage for PMN J2219-2719 is 
Ri2-9 keV. 

21 For the new Chandra observation of PMN J2219-2719 the 
total number of counts in a circular aperture extending from 1" 
to 4" is 140^2! an d the number of counts predicted by our PSF 
simulation is 102. 




Fig. 4. — Maximu m-likelihood reconstruction (see §2.1 
of ITownsIev eTaTI [2006) of the full band X-ray image of 
PMN J2219-2719. The colors yellow- red- blue indicate the decreas- 
ing order of photon count densities. Additionally a super-resolved 
VLA contour plot at 4.9 GHz is overlaid. The image resolution of 
the radio contour plot is 0.6" X 0.6". The restoring beam is the 
blue circle plotted at the lower-right corner. The contours begin 
at 0.2 mjy per beam and increase by factors of 2. 

tending from 1" to 4", which implies that the resolved 
X-ray emission has a count rate in the full band of 
(2.1 ± 0.7)% that from the core. Such a level of puta- 
tive X-ray jet emission is consistent with that estimated 
by Lopez et al. (2006) using the Chandra archival ob- 
servation of PMN J2219-2719. Additionally, our anal- 
yses suggest a south-western X-ray angular extension 
of 2" from PMN J2219-2719 (~ 14 kpc in the RLQ 
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rest frame), consistent with the possible ~l-2" exten- 
sion noted by Lopez et al. (2006). The total number of 
counts that lie at >1" from the core and that are inside 
the south-western semi-circular sector presente d in the 
upper-right inset of Figure [3] is 47^ (errors from lGehreld 
fl986h . In the same area the number of expected counts 
for a point source (from the MARX PSF) is 24. Based on 
a Poisson distribution with expected mean 24, the null 
hypothesis of obtaining 47 counts or more is 2.12 x 10~ 5 . 
Therefore the extension found in PMN J2219-2719 has 
a significance of « 4.1cr. The putative X-ray jet lumi- 
nosity in the 2-10 keV rest-frame of PMN J2219-2719 
is log L2-10 ~ 44.6 (<~2% of the total £2-10 luminosity 
of PMN J2219-2719; see §£2$. Given this high X-ray 
luminosity, the extended emission for PMN J2219— 2719 
is unlikely to have a non-AGN origin; e.g., starburst ac- 
tivity is typically fainter than log L2-10 ~ 42.5 (see, e.g., 
iBauer et al.ll2004l) . 

We find a low-level secondary radial structure in the 
4.9 GHz image, which can be seen in the south west- 
ern tip of the contour plots in the right panel of Fig- 
ure [2] This structure is apparent more clearly when 
we super-resolve the image using a restoring beam of 
0.6" x 0.6" (see the contour plot in Figure H). The 
south western extension at 4.9 GHz has a peak flux 
density of 0.54 mJy per beam, around 5 times higher 
than the rms of 0.11 mJy per beam. This 4.9 GHz 
extension in PMN J2219-2719 is spatially distinct 
from the southernly exten sion claimed at 1.4 GHz by 
ICheung et all (|2006l l2008[ ). At this southernly decli- 
nation, the VLA synthesized beam at 1.4 GHz has di- 
mensions of 2.70" x 1.00", with the major axis oriented 
north-south. Cheung et al. use a restoring beam of 
0.75" x 0.75" to deconvolve this source, super-resolving 
the major axis by almost a factor of four. The claimed 
secondary radio component extends to the south, in the 
direction of maximum super-resolution and highest un- 
certainty. We have similarly super- resolved the 1.4 GHz 
image of PMN J2219— 2719 and are unab l e to r e plicat e 
the component shown by ICheung et all (|2006l I2008D . 
To compare agreement between the X-ray extension 
and the secondary radio structure found at 4.9 GHz in 
PMN J2219— 2719, we perf ormed a maximum-lik elihood 
reconstruction (see §2.1 of iTownslev et al"1l2006l ) on the 
full band image (250 iterations). The reconstructed 
image shows various marginal features (1-3% of the 
core intensity) at all position angles, although somewhat 
stronger to the south and west. In Figure [4] we show a 
reconstructed image (yellow-red-blue indicate decreasing 
order of count density) with logarithmic intensity scal- 
ing. From Figure 0] we notice that the X-ray extended 
emission of PMN J2219— 2719 approximately matches 
the weak radio extension in the SW direction. 

Additionally, for the new Chandra observation of 
PMN J0235-1805, the excess of X-ray counts that lie at 
>1" from the core were consistent with zero when com- 
pared with the MARX PSF. The line-of-sight projected 
jet emission in PMN J0235-1805 and PMN J2219-2719 
is < 2.6% of the core emission in these sources. Simi- 
larly, previous work on the z = 4.30 RLQ GB 1508+5714 
dem onstrated that the jet-to-core X-ray flux ratio was ~ 
3% (|Siemiginowska et all 120031: lYuan et al.ll2003| ). The 
jet X-ray emission in PMN J2219— 2719 is also consistent 



with a mean of ~ 4% (a w 6%) jet-to-core X-ray flux ra- 
tio fojono^n_j 1 _sjmpkj[20 object of z < 2.6 quasar jets 
bv iMarshall et al.l (|2005l . 1201 111 . Therefore in neither of 
our new Chandra observations do we detect a fractionally 
large jet emission when compared to the core emission. 
As a consequence, we do not find any evidence of en- 
hanced jet emission at the high redshifts of these sources 
as we might ex pect from the IC/CMB model (see, e.g., 
ISchwartzll2002l ). 

3.2. Spectral analysis 

The spectral analysis was performed using xspec 
v.12.0; the fitting bands used were 0.5-8 keV for the 
Chandra observations, 0.5-10 keV for the XMM-Newton 
MOS camera observations, and 0.3-10 keV for the 
XMM-Newton pn camera observations. For Chandra 
we performed fits for each of the new and archival ob- 
servations. Since the sources PMN J0235— 1805 and 
PMN J2219-2719 have both new and archival obser- 
vations with Chandra, we have performed a joint spec- 
tral fitting of each new observation with its respective 
archival observation. For these joint fits all the spectral 
parameters excluding the normalization are constrained 
to be the same. For the XMM-Newton data we per- 
formed fits on each of the EPIC camera spectra sepa- 
rately; we also performed joint fits using the three EPIC 
cameras. The fits for each EPIC camera separately are 
consistent with the joint fits, and therefore we present 
only the joint fits in this work. 

One of the objectives of the spectral analysis is to 
determine if intrinsic absorption is important in typi- 
cal high-redshift RLQs. As a first approach to assess 
the presence of absorption in our sources, we group the 
spectra to have between 10-30 counts per bin in each 
observation/instrument. Using these grouped data we 
fit the observations (using the x 2 statistic) above 5 keV 
in the rest frame (about 1 keV in the observed frame) 
with a power-law and then extrapolate it to lower ener- 
gies. Through the use of this method, the residuals found 
at low energies may indicate absorption in the spectra. 
Evidence for absorption is only seen in the spectrum of 
SDSS J0011+1446 (see Figures [5] and E|) . The qualita- 
tive analysis above is also consistent with our spectral 
fits. These spectral fits were performed in the f ull fit- 
ting b and of each instrument using the C-statistic (|Cashl 
[1971 . 22 We first try a fit using a power-law model with 
Galactic absorption and intrinsic absorption (the xspec 
model PHABS*ZPHABS*POW). The Galactic absorption 
colum n densities were obtained from lDickev fc Lockmanl 
(1990); the errors on the Galactic column densities 
range between (1 - 3) x 10 19 cm" 2 (|Elvis et al.lll994ri 
Murp hy et al.fi .996). In all the fits solar abundances were 
assumed for the Galactic and intrinsic absorption. We 
do not find evidence of significant intrinsic absorption in 
five out of our six new RLQ observations. This can be 
inferred from the best-fitted parameters obtained from 
the absorbed power-law model. The column densities for 
these five sources are consistent with zero. Additionally, 
contour plots of T vs. N-r indicate that a null intrinsic 
column density is inside the 68% confidence region (see 
Figures [5] and O. The source that shows evidence for ab- 
sorption at >99% confidence is SDSS J0011+1446. The 

22 No spectral grouping was used for the C-statistic fits. 



TABLE 4 

Best-fit X-ray spectral parameters 





GALACTIC Nn a 7V H b 








OBJECT NAME 


OBSERVATION ID (10 20 cm -2 ) (10 22 cm -2 ) 


r 


/0.5-2 C 


log L 2 -io d C-stat(DOF) 



XMM-Newton observations 



SDSS J0011+1446 


0600090101 


4.3 


i 7-,+O.bV 

L - ' i -0.64 


1.87±0.07 


13.74±0.75 


46.62±0.02 


1146.4(1336) 


PMN J0214-0518 


0600090401 


2.6 


< 1.05 


2.10±0.10 


1.88 ±0.19 


45.55±0.04 


887.8(903) 


SDSS J0839+5112 


0301340101 


3.3 


< 2.33 


1.61±0.12 


4 42+ - 73 


45.99±0.06 


407.3(446) 


SDSS J1309+5733 


0301340501 


1.3 


< 2.21 


1 69+° 15 


1.68 ± 0.24 


45.53±0.06 


478.0(542) 



Chandra observations 



SDSS J0011+1446 


3957 


4.3 


< 6.1 


1.90±0.26 


12.31+^ 


46.58±0.10 


109.6(90) 


PMN J0214-0518 


4767 


2.6 


< 4.0 


n ic+0.53 
z - 1D -0.48 


n 70+1.23 

°- ' -0.87 


45.85±0.12 


24.1(37) 


PMN J0235-1805 


10306 


2.7 


< 4.3 


1.86±0.25 


2.32 ± 0.37 


45.71±0.07 


85.4(104) 


PMN J0235-1805 


4766 


2.7 


< 15.7 


1.70 


1 co+0.72 
L - OZ -0A6 


45.52±0.16 


23.3(18) 


PMN J0235-1805 


10306 and 4766 


2.7 


< 4.8 


1.84±0.24 




108.7(122) 


SDSS J0839+5112 


3562 


3.3 


< 3.4 


1.76±0.32 


5.04±1.04 


46.06±0.09 


54.1(71) 


SDSS J1309+5733 


3564 


1.3 


< 24.7 


1.70 


2.43±0.71 


45.71±0.13 


24.0(26) 


PMN J2219-2719 


10305 


1.4 


< 0.6 


1.34±0.06 




46.28±0.02 


443.9(508) 


PMN J2219-2719 


4769 


1.4 


< 3.2 


1.33±0.18 


' ■ io -0.94 


46.03±0.06 


154.2(159) 


PMN J2219-2719 


10305 and 4769 


1.4 


< 0.6 


1.34±0.06 




598.1(668) 



Note. — The best-fit column density, were obtained from a model consisting of a Galactic-absorbed power-law with intrinsic 
absorption at the redshift of the source; the same model has been used to estimate the best-fit photon index, and C-statistic for the 
observations of SDSS J0011+1446. The best-fit photon index, and C-statistic of all the other observations were obtained from a model 
consisting of a Galactic absorbed power-la w. Errors and upper limits represent 90% confidence intervals for each value, taking one 
parameter to be of interest (AC=2.71; e.g., |Avni||1976D. 

a Galactic absorption column densities (from Dickey &c Lockman 1990), obtained at the optical coordinates of the sources (see Tablcs[T] 
and [2j through the use of the HEASARC Nu tool at http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
k Intrinsic column density in units of 10 22 cm~ 2 . 

c Absorption-corrected flux in the observed 0.5—2 keV band in units of 10 — 14 erg cm~ 2 s _1 . This flux is obtained from the best 
fitted parameters of an absorbed power-law fit. Note that the intrinsic absorption is only important for SDSS J0011 + 1446; for the 
XMM-Newton observation of this source /o.5_2 — (14.16±0.75) X 10 -14 erg cm~ 2 s~ 1 when we do not correct for intrinsic absorption. 
^ Logarithm of the absorption-corrected luminosity in the 2 — 10 keV rest-frame band; the luminosity is in units of erg s — 1 . 



absorption found in this source is of moderate strength 
(TVh ~ 1-7 x 10 22 cm -2 ), in agreement with the absorp- 
tion found in other high -redshift RLQs (e.g.. iPage et al.l 
120051: lYuan et al.ll2006l) . To verify that the absorption 
detected cannot be reproduced by unaccounted uncer- 
tainties in the Galactic absorption, we have fitted the 
XMM-Newton spectra (joint fit of the EPIC cameras) of 
SDSS J0011+1446 using a Galactic absorbed power-law 
model (i.e., the xspec model phabs*POw). The best- 
fitted absorption obtained, with the absorption at z = 
as a free parameter, is (7.7 ± 1.5) x 10 20 cm -2 (error 
at the 90% confidence level), well above the expected 
value for Galactic absorption (see Table 2]) . Based on 
the fact that just one observation presents absorption, 
we re-fit all the other observations (new and archival) 
with a power-law model with only Galactic absorption 
(the xspec model PHABS*POW). The best-fit spectral 
parameters (errors and upper limits at the 90% confi- 
dence level) are presented in Table [4] In this table, the 
best-fit photon index and C-statistic, of all the observa- 
tions that do not present absorption, are obtained from 
a model consisting of a Galactic-absorbed power-law. In 
the case of SDSS J0011+1446 the best-fit photon index 
and C-statistic presented in Tableware obtained from a 
model consisting of a Galactic-absorbed power-law with 
intrinsic absorption at the redshift of the source. The 
mean photon index for our sample is 1.74 ±0.11 (la error 
on the mean) , in agreement wi t h samples of low-redshif t 
RLQs (jSambruna et al.l 119991: iReeves fe Turner! |2000D . 
Based on this, in Table U for the Chandra archival ob- 



servations of PMN J0235-1805 and SDSS J1309+5733 
we have fixed the power-law photon index to T = 1.7 in 
order to avoid degeneracies associated with the limited 
numbers of counts (< 30; see Table [5]) in these observa- 
tions. 

Since the average redshift of our sources is close to four, 
the rest-frame band of the fitting for each source reaches 
energies up to ~ 50 keV. Therefore, the source spectra 
are sampling any potential Compton-reflection compo- 
nent which peaks at rest-frame energies of ~ 30 keV (e.g . , 
iLightman fe White! [ilSM iMagdziarz fe Zdziarskilll995l) . 
In general we would not expect reflection signatures to 
be strong in these sources due to likely dilution by the 
jet-linked X-ray continuum. To check for reflection we 
fitted our spectra above 5 keV in the rest-frame using 
the PEXRAV model. From these fits, we find that the 
reflection parameter is consistent with zero in all our ob- 
servations, and therefore we do not find that reflection is 
significant in any of our sources. We also used Gaussian 
line profiles to fit our spectra in the energy range where 
the Fe Ka line is found (6.4 keV rest-frame for neutral 
iron). Through this procedure we find that the Fe Ka 
best-fitted equivalent widths (EWs) are consistent with 
zero for all our observations, and therefore there is no ev- 
idence of Fe Ka signatures in our spectra. For the new 
observation of SDSS J0011+1446, which is the observa- 
tion of highest S/N, we obtain that the EW of Fe Ka 
is < 170 eV. The upper limit on the EW of Fe Ka line 
for the other sources is less well constrained. The ab- 
sence of significant Fe Ka emission lines is consistent with 
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Fig. 5. — X-ray spectral fits of a power-law with Galactic absorption to the Chandra observations of PMN J0235— 1805 and 
PMN J2219— 2719. The best-fit power-law to the data above 5 keV rest-frame and extrapolated downward (solid line) is shown. The 
open diamonds correspond to the data points of our most recent Chandra observations of PMN J0235— 1805 and PMN J2219— 2719. 
The open squares correspond to the data points of the Chandra observations of PMN J0235— 1805 and PMN J2219— 2719 analyzed in 
I Lopez ct al. (2006). The Ax residuals show the deviation of the data from the model in units of a with error bars of size unity. In each 
residual panel we have marked with a dotted vertical line the 5 keV rest-frame energy. The inset in each panel shows 68%, 90%, and 99% 
confidence contours for the intrinsic absorption (Afjj) and photon index (r) using the simultaneous C-statistic fits of observations 4766 and 
10306 for PMN J0235-1805 and 4769 and 10305 for PMN J2219-2719. 



TABLE 5 

TWO EPOCH X-RAY FLUX VARIABILITY TEST. 





/0.5-2 a 


/().5-2 a 


At h 






OBJECT NAME 


(epoch 1) 


(epoch 2) 


(days) 


x 2 


%sig c 



SDSS J0011+1446 


12 S1+ 2 - 04 


13.74±0.42 


406 


0.5 


51.1 


PMN J0214-0518 


o 79 +0.77 
J - ' -0.52 


1.88±0.11 


446 


11.98 


>99.9 


PMN J0235-1805 


, „+0.46 
1 - OZ -0.27 


2.32±0.24 


347 


2.4 


87.8 


SDSS J0839+5112 


5.04±-0.64 


4 42+ - 47 
4 - 4z -0.32 


218 


0.6 


56.8 


SDSS J1309+5733 


o 4^+0.39 


1.68±0.16 


193 


2.3 


87.0 


PMN J2219-2719 


7.15±0.59 


12.78t°il 


459 


54.8 


>99.9 



a Absorption-corrected flux in the observed 0.5 — 2 keV band in units of 
10 -14 erg cm -2 s _1 . For each source the fluxes at epoch 1 and epoch 2 cor- 
respond to those obtained from the Chandra archival observations and from 
the new observations (either with Chandra or XMM-Newton), respectively. 

Rest-frame timcscalc. 
c Percentile of significance. 



what is expected f or luminous (Lx > 10 46 erg s Ta- 
ble H I) quasars (e.g., iNandra et atlll997tlPage et al.ll2004t 
iBianchi et alll2007f l. 

3.3. Short and long-term variability 

We searched for rapid variability in the new 
XMM-Newton and Chandra observations by applying 
the Kolmogorov-Smirnov test to the photon arrival 
times; however, no significant variability was found. 



Long-term variability was assessed by comparing the 
fluxes in the 0.5-2 keV band from the archival obser- 
vations with those of the new observations. Comparing 
the fluxes in the two observation epochs of each source 
(Table |4|) , we find that one third of our sample (2 out 
of 6) shows potential variability. The sources which 
possibly present variability are PMN J0214— 0518, and 
PMN J2219-2719. In order to test the significance of 
this variability, we calculate for each source the \ 2 statis- 
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Fig. 6. — X-ray spectral fits of a power-law with Galactic absorption to our new XMM-Newton observations of z > 4 RLQs. The 
best-fit power-law to the data above 5 keV rest-frame and extrapolated downward in energy (solid line) is shown. The Ax residuals show 
the deviation of the data from the model in units of a with error bars of size unity. The open diamonds, open squares, and filled squares 
correspond to the data points of the PN, MOS1 and MOS2 spectra, respectively, for our XMM-Newton observations. In each residual 
panel we have marked with a dotted vertical line the 5 keV rest-frame energy. Note the systematic negative residuals below 0.8 keV 
for SDSS J0011+1446, corresponding to the detected intrinsic absorption. The inset in each panel shows 68%, 90%, and 99% confidence 
contours for the intrinsic absorption (Nn) and photon index (T) using the simultaneous C-statistic fits of the pn, MOS1 and MOS2 data 
points. 



tic where the data points are the fluxes of each observa- 
tion (with their respective la errors) and the model is a 
constant flux with the best fit of the two epochs. The 
X 2 value provides a statistical test of the null hypothesis 
that the flux of each epoch is equal to the best-fitted flux 
of the two epochs. In Table [5] we list the fluxes for each 
source in the Chandra archival observations (epoch 1) 
and in the new observations (epoch 2). In this table we 
also have listed the values of \ f° r our model to test vari- 
ability. For PMN J0214-0518, and PMN J2219-2719 we 



reject the null hypothesis with > 99% significance. 23 The 
rest-frame time scales for the variability are 446 days and 
459 days for PMN J0214-0518 and PMN J2219-2719, 
respectively. The amplitude factors of the flux vari- 
ability are «2.0 and for PMN J0214-0518 and 
PMN J2219— 2719, respectively; these are in agreement 
with the typical amplitude fa ctors of lower redshift RLQs 
(e.g JBrinkmann et alJll997f ). The two sources that show 

23 We arrive to the same conclusion if we compare the fluxes in 
the 0.5—8 keV band from the archival observations with those of 
the new observations. 
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Fig. 7. — Radio-loudness parameter versus best fitted column density (Nn; left panel) and photo n index (T ; right panel) for our 
sample (stars) and z > 2.0 RLQ observation s with Chandra and XMM-Newt on (circles a re from Page et al.j 2005, triangles from 
IShemmer et al.l 2005, squares from lYuan et all 2006, and inverted triangles from Young et al. 2009). The arrows in the lef t panel cor - 
respond to column density upper limits. The dotted line is obtained by using the the IDL Astronomy library tool linmix_err l lKellv|[2007l ) 
on our enlarged sample. 



evidence for variability are the two with the longest time 
scales between observations, but this may be merely co- 
incidental. 

3.4. Power-law and absorption correlation tests 

In order to investigate possible relations between X-ray 
spectral parameters, redshift, and Rl, we have enlarged 
our original sample using some of the latest high-redshift 
(z > 2) RLQ X-ray measurements. The Rl vs. z dia- 
gram of the enlarged sample can be seen in Figure [JJ 
in which the sources from each survey have been indi- 
cated with different symbols. We have also plotted in 
Figure [7] diagrams of Rl vs Nh (left panel) and Rl vs. 
r (right panel) for RLQs at z > 2. From the right panel 
we infer a likely decrease of T with increasing Rl- To 
quantify possible dependences between the parameters, 
we have calculated Spearman correlation coefficients of 
T vs. Rl, r vs. z, Nh vs. Rl, and also ATh vs. z; 
the results are given in Table |5] considering sources with 
z > 2 and z > 3 separately. We find an anticorrela- 
tion between T and Rl ; this anticorrelation is significant 
(at the 99.9% confidence level) for RLQs with z > 2 and 
with z > 3. The hardening of the spectra with increasing 
Rl indicates a rising jet contribution to the X-ray spec- 
tra as a function of Rl- This conclusion is reinforced 
by the fact that the excess X-ray luminosity in RLQs 
(with re spect to RQQs) inc reases with Rl as quantified 
by, e.g., iMiller et all (|2011[ ). Notice that, with the ex- 
ception of PMN J2219— 2719, our sources possess mod- 
erate values of Rl (see Table [3]). The photon index for 
PMN J2219-2719 of T « 1.3 is clearly harder when com- 
pared with the photon indices of the other RLQs in our 



sample (on average » 1.7; see 33. 2j) : this is in agreement 
with the trend found. Additionally, the soft spectral in- 
dex (r « 2.1) of PMN J0214-0518 (log R L ~ 2) is also 
in accordance with the trend found (see Figure [7]) . 

The Pearson correlation coefficient of T with log Rl for 
our extended sample is —0.487 (significant at the 99.7 % 
confidence level); therefore, as a second step we fit the re- 
lation between T and log Rl with a simple linear model. 
This fit is not statistically acceptable because the error 
bars on T are not representative of the intrinsic disper- 
sion of the data around a possible T — log Rl relation. In 
view of this, we adopt th e IDL Astronomy Library tool 
LINMIX_ERR (|Kellvll2007T ) to estimate the parameters de- 
scribing a linear model between T and log Rl (i.e., find 
a and /3 of T = a + f3 x log Rl)- This IDL routine, which 
accounts for measurement errors, nondetections, and in- 
trinsic scatter, is implemented with a Baycsian approach 
to compute the posterior probability distribution from 
the observed data. From the random draws (2 x 10 5 re- 
alizations) from the posterior probability distribution we 
find that the median values of the parameters a and (3 
are given by 

T = (2.03 ± 0.11) + (-0.15 ± 0.04) x log R L . (1) 

The r — log Rl relation has alre ady been found in the 
past in z 5, 2 samp les of RLQs (e.g.. lWilkes fc Elvislll987l : 
iReeves fc Turnerl uOOO) . To test if the relation found in 
equation [JJ is consistent with that fo und at lower red- 
shifts, we selected z < 2 RLQs from IReeves &: Turner! 
( 2000). Using linmix_err (2 x 10 5 realizations) we find 
that the median values of the parameters a and f3 are 
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Fig. 8. — Contour plots of a vs. for the model V = 
a + ft X log; Rl usin g the the IDL Astronomy Library tool LIN- 
mix_err QKellvll200"7l ). The solid and dotted contours are the 68% 
and 99% confidence regions of 2 X 10 realizations for our and the 
Reeves & Turner data, respectively. 



TABLE 6 

Correlation analysis results for a combined sample of 
high-redshift radio-loud quasars. 



Tested parameters Redshift range 7V a 



rs 



%sig 



r vs. r l 
r vs. r l 
r vs. z 
r vs. z 

N H vs. R L d 
N n vs. R L d 
7V H vs. z d 
TVh vs. z d 



> 2.0 

> 3.0 

> 2.0 

> 3.0 

> 2.0 

> 3.0 

> 2.0 

> 3.0 



35 
22 
35 
22 
35 
22 
35 
22 



-0.542 
-0.637 
0.063 
0.341 
-0.165 
-0.344 
0.370 
0.395 



99.9 
99.9 
28.1 
88.0 
66.4 
88.5 
97.0 
93.0 



Note. 



The sample used to calculate the correlation parame- 



ters corresponds to radio -loud quasa r s presented in Tab les [Tl and [2] 
and those extracte d from Page et al.l 120051. Shc mmer et all (2005), 
lYuan et all 120061) . and lYoung et alj l^OOsh . 
a Number of RL quasars in each redshift bin. 

Spearman correlation coefficient. 
c Percentile significance of the correlation. 

d The correlation coefficients involving TVh were calculated using 
survival analysis with the ASURV code (Isobc et al. 1986). 

given by 

T = (2.09 ± 0.16) + (-0.12 ± 0.05) x log R L . (2) 

The linear para meters obtained for the z < 2 sources of 
I Reeves fc Turner! (equation [2]) are consistent with those 
found for our sources (see equation 1, and Figure [8]); 
therefore, we do not find evidence that V— log Rl relation 
changes with redshift. 

The correlation tests in Table [6] give indications that 
Nh could be in creasing with z for sources with z > 2 (i n 
agreement with lReeves fc Turnerl2000HPage et"aT]|2005D . 
The correlation of Njj with z is at the ~ 97% confidence 
level. Therefore, an apparent correlation between A^h 
and z may be the result of observational bias, since the 
threshold of absorption that can be measured increases 
with redshift. 24 

24 The minimum rest-frame energy of the observed X-ray band 



4. SUMMARY AND CONCLUSIONS 

We have reported the X-ray properties of six high- 
redshift (z > 4) RLQs, observed with XMM-Newton 
(four observations) and Chandra (two observations) as 
a follow-up of earlier Chandra snapshot observations. A 
key aspect of our observed sample is that it spans lower 
values of radio loudness than past studies of high- redshift 
RLQs in the X-ray regime; the moderate radio-loudness 
objects studied are much more representative of the over- 
all RLQ population. The main conclusions from our 
spectral, imaging, and variability analyses are the fol- 
lowing: 

1. The power-law X-ray continua of our targets have 
a mean photon index of (r) = 1.74 ± 0.11, con- 
sistent with measurements of lower redshift RLQs. 
They also follow an anti-correlation between pho- 
ton index and radio loudness seen at lower red- 
shifts. Using linear fits to the T — log Rl relation 
for an extended sample at z > 2, we find consis- 
tency of this dependence in relation to lower red- 
shift (z < 2) samples. We find no evidence of 
Compton-reflection continua or iron Ko emission 
lines. The measured continuum shapes, combined 
with the fact that our targets are on average « 6 
times brighter in X-rays than RQQs of similar UV 
luminosity, indicate that jet-linked emission domi- 
nates the X-ray continua. 

2. We find evidence that one of our targets, 
SDSS J0011+1446, has significant X-ray absorp- 
tion with a column density of N-r ~ 1.7 x 
10 22 cm~ 2 . We set useful limits upon X-ray ab- 
sorption for the rest. The incidence of X-ray ab- 
sorption in our sample appears roughly consistent 
with that for more highly radio-loud objects at sim- 
ilar redshifts. 

3. In the Chandra observation of PMN J2219-2719 
we detect a likely X-ray jet with an extension of 
~ 14 kpc. The X-ray luminosity of this putative jet 
is w 2% that of the core. In the Chandra observa- 
tions of PMN J0235-1805 and PMN J2219-2719 
we do not find evidence of enhanced extended 
X-ray jet emission as we might expect from the 
IC/CMB model. We do not find any evidence of ra- 
dio extension in the 1.4 and 4.9 GHz VLA images of 
PMN J0235-1805. We find a secondary structure 
in the 4.9 GHz VLA image of PMN J2219-2719 
that matches the X-ray extension found in this 
source; due to low resolution we are unable to see 
this structure in the 1.4 GHz VLA image. 

4. Two of our targets, PMN J0214-0518 and 
PMN J2219-2719, show evidence for X-ray vari- 
ability by 80-100% on rest-frame timescales of 
450-460 days. 

increases as (1+z), e.g., for z ~ 4 the minimum rest-frame energy 
of the observed band is > 2 keV. A typical intrinsic absorption of 
column density ~ 10 22 cm -2 shows a signature in the spectra at 
rest-frame energies < 2 keV. Therefore, the threshold absorption 
that will produce a significant feature in the spectrum increases 
with redshift. 
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